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^ . Abstract 

We consider the non-radiative two body decay of a neutrino to a daughter 
neutrino with degraded energy and a very light particle (Majoron). Ultrahigh 
energy neutrinos from an astrophysical source like a Gamma Ray Burst under- 
going this decay process are found to produce different number of events in the 
■ detector depending on whether they are Majorana or Dirac particles. The next 

generation large scale experiments like Multi-OWL is expected to provide us 
an accurate determination of the flux of neutrinos from astrophysical sources 
00 ' and this may enable us to distinguish between the Dirac and Majorana nature 

of neutrino. 
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The recent obsevations from Super-Kamiokande [I], KamLAND [2] and Sudbury 
Neutrino Observatory (SNO) [Sj have established mixing between different neutrino 
flavours. There are, however, several questions which remain unanswered. The con- 
ventional model leads to three flavour states arising out of three mass eigenstates 
with Am] 2 k 7 x \Q~ 5 eV 2 and Am| 3 k 3 x 10~ 3 eV 2 . The Majorana/Dirac nature 
of the neutrino and the stability of all the states are still open questions. The Majo- 
rana/Dirac nature of neutrino can be resolved unambiguously if neutrino less double 
beta decay is observed. This difficult experiment is being carried out at various lab- 
oratories. Other alternative methods to establish Dirac/Majorana nature of neutrino 
would be desirable and this work is an attempt in that direction. 

We have been motivated by the recent work of Beacom et al. [HE], who have 
considered the possibility that some of the neutrinos could be unstable. Depending 
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upon the decay pattern they have noted that high energy neutrinos coming from 
astrophysical sources can have different flavour ratios. Such a decay of a "heavy 
neutrino" to a "light neutrino" could occur say by an emission of a spinless particle. 
For example, if neutrino mass pattern is given by normal hierarchy and if only the 
lightest neutrino is stable then the flavour ratio of the astrophysical neutrinos reaching 
the Earth would be v e : : v T — 6 : 1 : 1. An experiment like IceCube which 
could distinguish different neutrino flavours, may provide observational evidence for 
the different decay scenarios. 

One has to consider detecting very high energy (> 10 15 e^) neutrinos to ensure the 
background from atmospheric neutrinos are negligible. Astrophysical sources which 
could emit high energy neutrinos are Active Galactic Nuclei (AGN) and Gamma Ray 
Bursts (GRBs). Neutrino fluxes from such sources have been studied considering 
different theoretical models, and are expected to produce large number of events 
in experiments being planned e.g. Antarctic Impulsive Transient Antenna (ANITA) 
[7j, Orbiting Wide-angle Light-collector (OWL) 8j. In this work, we explore the 
possibility of distinguishing the Dirac/ Major ana nature of the neutrinos from these 
observations. The basic idea is this : If a "heavy neutrino" decays to a "light neutrino" 
with a spin flip then the light neutrino will be sterile if the neutrino is a Dirac particle. 
However, if the daughter neutrino is a Majorana particle, this would still be an active 
neutrino and would produce a signal in the detector. We consider the neutrinos from 
a single Gamma Ray Burst. The nature of Gamma Ray Bursts, the emitted photon 
spectrum in the keV-MeV energy range and the subsequent emission of afterglow 
have been subjects of intense study in the last few years. Based on a "fireball" model 
of GRB, the spectrum of high energy neutrinos from all GRBs has been obtained in 
The expected neutrino spectrum from a single GRB in the same model has been 
calculated in [TTH ITT]. 

We show that the number of active neutrino events would be different in the de- 
tector depending on whether neutrinos are Dirac or Majorana particles. In particular 
we have looked at upward going tau neutrinos having energies above 3 x 10 15 ey tak- 
ing into account of their regeneration as they traverse the Earth and look for signals 
of Dirac vs Majorana nature. OWL [S] experiment will use the Earth's atmosphere 
as a detector. Multi-OWL 12 would cover the Earth's horizon with a field of view 
of 120° as it would view from an orbit of about 1000 Km. The area on the Earth's 
surface, which would be within the field of view of Multi-OWL, would be used as 
the effective target area for the upward going neutrinos. Further, this experiment is 
suited for detecting tau neutrinos through observation of ultraviolet (UV) fluorescent 
light emitted by secondary charged particles. For energies above 10 15 eV the main 
contribution would be from v T for the upward vertical showers as v e and z/ M would 
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be absorbed in matter. The prospects of neutrino detection above energy 10 15 eV^ by 
OWL experiment have been discussed in |12j . 

In order to obtain results which do not depend sensitively on the model dependent 
parameters of the GRB, we have attempted to consider ratios of events belonging 
to different energy ranges. We discuss about the experimental accuracy needed to 
understand whether neutrino is a Dirac or Majorana particle based on the event rates 
calculated for Multi-OWL experiment. 



2 The flux of neutrinos from a single Gamma Ray 
Burst 

We consider a typical Gamma Ray Burst (GRB) whose parameters are inferred from 
observations, as our source of UHE (ultrahigh energy) neutrinos and calculate the 
number of tau neutrino events from it in Multi-OWL. There are discussions on the 
parameters of a GRB and the mechanism of neutrino production in GRBs in the 
fireball model in APPENDIX A. We have used the neutrino spectrum from a single 
GRB in the fireball model |TT1 I13j. The ultrahigh energy neutrino spectrum in the 
source (GRB) rest frame is as given below. 

dN (E) if ( E "/ E t) al E » < E l 

^^oc-M (E u /Eir El<E v <K (1) 

dK K { {E v /ElY\E v /ElY 2 E v > El 

where, ll E h u n is the neutrino spectrum break energy, "al" and "a2" are the neutrino 
spectral indices below and above the neutrino spectrum break energy {E^) respec- 
tively, is due to the break energy in the low energy photon spectrum of a GRB 
(see APPENDIX A for discussion). u E^ n is the break energy in the neutrino spectrum 
due to the synchrotron energy loss by high energy pions. 

The final neutrino spectrum on Earth can be derived in the following way. As- 
suming that the GRB is occuring at a redshift "z" , the energy per unit time and unit 
area to be measured by a detector on Earth is 



Fuohs ~ 4ird 2 L (2) 

where, "g?l" is the luminosity distance of the source. The luminosity in the source 
rest frame ll L v n can be expressed in terms of the number of neutrinos per unit energy 
and unit time in that frame. 

L ^f- E ,^A dK (3) 
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where, dN ^ v ^ = ^ d ^J^df"^ dt v , "£</' is the duration of the burst in the source rest 
frame. We denote the observed neutrino spectrum on Earth per unit area, unit energy 
and unit time by dM d ^ u ° h ^f u ° ba ^ • From eq.(2) and (3) we can write 

P dM' vobs {E vobs ,t uobs ) ^ 1 r^W dK(Ej 

V ° bs IF ~M dtj vobs = -7-J2 / "~Tp~~JT~ v \ ' 

Hence using E uobs = E v j (1 + z) and di = (1 + z)d z , 

dM' uobs {E vobs ,t vobs ) 1 dN' u {E U) t u ) 



dE uobs dt uobs Aa[d\ dE v dt v 



(5) 



u d z " is the comoving distance. The number of neutrinos received on Earth per unit 
area and unit energy is 

dM uobs (E uobs ) _ 1 , dN u (E v ) 

dE vobs ~ 4irdl [ +Z) dE v [ } 



where, dt vobs = (1 + z)dt u , ^-f = \t "°dt dt vobs and t dobs = t d (l + z). 

In a spatially flat universe we calculate the comoving distance of a source using 
— 0.73, Q m = 0.27 and H Q = llkmsec -1 Mpc~ x from [T5]. The source spectrum 
contains v e and in the ratio 1:2. As the electron and muon neutrinos propagate 
from the source to the observer, the electron neutrinos oscillate to muon neutrinos 
which in turn can oscillate to tau neutrinos. If we assume U e s = 0, the solar neutrino 
mixing angle 8 & to be 30° and the atmospheric neutrino mixing angle 8 atm to be 45°, 
then due to the effect of neutrino oscillations alone, we expect an equal number of 
z/ e , and v T on Earth. Inclusion of the neutrino decay scenario causes the flavor 
ratio of the three active neutrino species to deviate from 1:1:1 4j. We assume the 
decay scenario of neutrinos where one neutrino decays to a daughter neutrino and 
a very light or massless particle. The daughter neutrinos do not disappear and the 
neutrino mass spectrum is hierarchical. The energy dependent neutrino decay term 
is proportional to exp(—dm u /E u t), where u d" is the distance of the source from the 
detector, u m u " is mass of neutrino and "t" is its rest-frame lifetime. We assume the 
neutrino decays to be complete when we receive them on Earth. This assumption is a 
reasonable one because there is a lower limit on the value of djE v set by the shortest 
distances of the sources (typically hundreds of Mpc) and the maximum energy that 
will be visible in OWL. The daughter neutrinos would be slightly deflected from 
the direction of propagation of the parent neutrinos. Since, we assume that the 
source is emitting isotropically, the flux recieved on Earth would not change due to 
the angular deflection of the decay products. If neutrinos are Dirac particles then 
daughters, which flip helicity, are sterile and undetectable. In the case of Majorana 
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neutrinos, neutrinos and antineutrinos of both the helicity states are active and there 
is no distinction between a neutrino and an antineutrino. The expression for the 
neutrino spectrum after the appearance of daughter neutrinos with degraded energies 
is given in eqn.(7) of j3j. The energy dependence of the decay rate of neutrinos has 
been explicitly given in Appendix B. W(E', E) is the energy distribution function of 
the daughter neutrinos where, E' and E are the energies of the parent and daughter 
neutrinos respectively. This distribution function is different for Dirac and Majorana 
nature of neutrinos because in the Dirac case the neutrinos with flipped helicities 
are sterile in contrast to the Majorana neutrinos. We have also verified that in the 
case of radiative decay of neutrinos the energy dependence of the decay rate remains 
identical to the scalar decay scenario of neutrinos if one considers a coupling of the 
type a^Ff". 

3 Propagation of UHE neutrinos through the 
Earth 

The propagation of v T through the Earth's matter at very high energy does not 
diminish their number [IB] . The vertically upward going u e s and u^s are absorbed 
in Earth's matter above 10 16 eV and can not produce signals in the detector. This is 
because the lifetime of /x lepton is large (10~ 6 sec) and it loses all its energy before 
decaying. This is in contrast to a r lepton whose lifetime is 10 _13 sec, hence it decays 
much before energy loss occurs due to electromagnetic interaction. The z/ r s produce 
r leptons in charged current interactions with matter nuclei. A r lepton can decay in 
many ways but in all of these decay processes another v T is produced in the final state. 
The v T regenerated in this way carries one third of the energy of the primary v T |17j . 
Also there are neutral current interactions and the v T generated in these interactions 
carries about half of the primary v T energy [TBj. High energy i> T s initiate a series of 
interactions in the Earth and ultimately there is generation of v T s of lower energies. 
Once the mean free path of interaction of u T with nucleon of Earth's matter becomes 
comparable to the radius of the Earth, it directly propagates to the detector. There 
are discussions on propagation of tau neutrinos through the Earth and their event 
rates in detectors in a series of recent papers [TBI EHl I2S ED 1221 • The procedure 
we have followed to calculate the probability of r lepton emission from the surface 
of the Earth has been discussed in APPENDIX C. The area on the surface of the 
Earth, which would be covered by the field of view of Multi-OWL, would act as the 
effective target area for the upward going tau neutrinos. One may look up ^2] for 
detailed discussion on how the larger field of view of Multi-OWL compared to OWL 
experiment would increase the effective target area on the Earth's surface for the 
upward going neutrinos. The r leptons emitted from the surface of the Earth would 
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produce UV fluorescent light in air while passing through the atmosphere. If they 
decay before reaching the atmosphere, there will be charged particles in their decay 
products which can produce UV fluorescent light in the atmosphere. UV fluorescent 
light will be emitted isotropically, hence OWL would be able to detect the charged 
particles from any direction [23]. The GRB v T spectrum on Earth folded with the 
probability of r emission from the surface of the Earth, calculated in APPENDIX C 
and the target area of the upward going neutrinos on the Earth's surface gives the 
number of r signals in the detector. 

4 Results 

The GRB parameters that have been used in our calculations are the following: 
a Lorentz factor T = 300, total energy emitted in neutrinos 10 51 erg, low en- 
ergy photon luminosity L 7 = 10 51 erg/sec, neutrino spectral indices al = 1.5 and 
a2 = —0.3,0.05,0.3 in eqn. (1), variability time t v = 10~ 3 sec, e e = 0.5, e# = 0.05. 
The values of these parameters have been inferred from low energy 7-ray observa- 
tions on GRBs. We have calculated neutrino events in Multi-OWL from a GRB at 
a redshift z = 0.5 and at a zenith angle 9 = 180°. Since, the neutrino spectrum 
break energy E h u is 2.26 x 10 15 eV^ for the GRB parameters used in our calculation 
and we are considering the neutrino events in the detector above energy 3 x 10 15 el / , 
the GRB neutrino spectrum above E h v only contributes to the number of events in 
the detector. The numbers of neutrino events calculated for Multi-OWL above en- 
ergy 3 x 10 15 eK are of the order of hundred. If we take the ratio of the numbers of 
events in the detector above different energies our results become independent of the 
parameters of the GRB e.g. T, z, e^ MeV and depend only on a 2 j2H]- We calculate 
the expected numbers of active neutrino events in the detector above energies 3 x 10 15 
and 2.2 x 10 16 eV^ for the two cases: (i) neutrino is a Dirac particle, (ii) neutrino is a 
Majorana particle. The ratios of the numbers of events above different energies are 
useful quantities for our study. They can be compared in future with their measured 
values in experiments like Multi-OWL. All GRB parameters except a 2 have been kept 
fixed for our results displayed in TABLE I, TABLE II, TABLE III and TABLE IV. 

TABLE I: Ratio of the numbers of vertically upward going v T events in Multi- 
OWL above energies 3 x 10 15 eV and 2.2 x 10 16 eV^ in the detector with mass ratio of 
daughter to parent neutrino 0.001 and a 2 = —0.3. 
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Unstable 


Branchings 


Ratio 


Ratio 

dM 




^2,^3 


irrevalent 


27.53 


29.99 


0.09 


^3 


#3-1 = 1 


24.92 


25.38 


0.02 


v-i 


#3->2 = 1 


25.84 


27.04 


0.05 


"3 


#3^1 = 0.5 

5 3 ^ 2 = °- 5 


25.4 


26.26 


0.03 



In TABLE I the unstable neutrino eigenstates/eigenstate and the decay modes con- 
sidered have been mentioned in first and second columns respectively. The third and 
fourth columns show the ratios of the numbers of neutrino events in the detector 
above energies 3 x 10 15 and 2.2 x 10 16 eU for Dirac and Majorana nature of neutrinos 
respectively. The fifth column displays by what fraction the ratios differ in the two 
cases. The ratios of the numbers of active neutrino events in the detector for the 
Dirac and Majorana nature of neutrinos differ by 9% for the decay scenario when 
both and v% are unstable and by 2% for -63^1 = 1. We have varied a 2 to study its 
effect on the ratios of the numbers of events above different energies and our results 
have been presented in the next two tables. 

TABLE II: Ratio of the numbers of vertically upward going v T events in Multi- 
OWL above energies 3 x 10 15 eU and 2.2 x 10 16 eU in the detector with mass ratio of 
daughter to parent neutrino 0.001 and a 2 = 0.05. 



Unstable 


Branchings 


Ratio 


Ratio 


{R™ - Ry T )/Ru T 


^2,^3 


irrevalent 


27.41 


29.85 


0.09 


^3 


£3-1 = 1 


24.83 


25.28 


0.02 


^3 


-83^2 = 1 


25.73 


26.93 


0.05 


^3 


5 3 ^i = 0.5 
B 3 ^ 2 = 0.5 


25.3 


26.15 


0.03 



TABLE II shows the ratios of the numbers of events in the detector above energies 
3 x 10 15 and 2.2 x 10 16 eU for a 2 = 0.05. 



TABLE III: Ratio of the numbers of vertically upward going v T events in Multi- 



7 



OWL above energies 3 x 10 15 eV A and 2.2 x 10 16 eV^ in the detector with mass ratio of 
daughter to parent neutrino 0.001 and a 2 = 0.3. 



Unstable 


Branchings 


Ratio 


Ratio 

dM 
K u T 


(R% ~ Ry T )/Ry T 




irrevalent 


27.31 


29.71 


0.09 




£3-1 = 1 


24.74 


25.19 


0.02 




-83^2 = 1 


25.63 


26.82 


0.05 


"3 


B 3 ^ = 0.5 
B^ 2 = 0.5 


25.2 


26.05 


0.03 



If we compare TABLE I, TABLE II and TABLE III, we find that the percentages 
by which the ratios of the numbers of events in the Dirac and Majorana case differ for 
different decay modes are independent of the value of ot 2 in the range —0.3 < a 2 < 0.3. 
Also, the ratios of the numbers of neutrino events do not vary significantly as we vary 
a 2 . We can possibly observe neutrino events from a large number of GRBs whose 
values of spectral index a 2 are in a range such that {R*f - B% )/B% are equal for 
all of them. In principle, if {R^f — Rf?) is greater than the sum of the theoretical 
and experimental errors in the ratios (R® or Rff) and we find that for most of the 
GRBs the experimentally measured ratios of the numbers of neutrino events above 
energies 3 x 10 15 and 2.2 x 10 16 eV A are closer to R® or Rff , then we would be able 
to say whether neutrino is a Dirac or a Majorana particle. We have discussed about 
the percentages of error in our theoretical calculation and neutrino experiments in 
section 5. 

We also study the effect of the mass ratio of daughter to parent neutrino on the 
number of events in the detector. We denote the masses of the daughter and parent 
neutrinos by mi and nth respectively. The area covered by Multi-OWL on the Earth's 
surface would be 3.3 x 10 13 m 2 [12]. Our results are displayed in TABLE IV. 

TABLE IV: Number of vertically upward going v T events in Multi-OWL from 
neutrinos having energies above 3 x 10 15 e\^, v 3 unstable, -63^1 = 0.5, -83^2 = 0.5 and 
a 2 = -0.3,0.3. 
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mi/m h 


Dirac 


Majorana 


Dirac 


Majorana 




«2 = —0.3 


«2 = —0.3 


a2 = 0.3 


«2 = 0.3 


0.01 


698 


735 


891 


938 


0.1 


699 


736 


892 


940 


0.25 


702 


744 


896 


949 


0.5 


711 


773 


908 


987 


0.75 


731 


829 


933 


1058 



The number of events increases with increasing neutrino mass ratios and this increase 
is more pronounced in case of Majorana neutrinos. 

5 Error Analysis 

The values of the parameters of a GRB are determined from 7-ray and afterglow 
observations on that GRB. These values inferred from observations certainly carry 
some errors due to the errors in the measured data by the experiments. They would 
be used in the theoretical prediction of the number of neutrino events from that 
GRB and the errors in the measurement of the values of the GRB parameters would 
propagate to the theoretical calculation of expected number of neutrino events. We 
have calculated the ratios of the numbers of neutrino events above different energies 
in the detector, which depend only on the value of the GRB parameter a^- If it is 
possible to determine «2 accurately from the low energy 7-ray spectrum emitted by 
the GRB then the error in our calculation would become negligible. For example, 
if the error in the measurement of «2, Aa 2 is 0.005 or 0.001 then the theoretical 
values of the ratios of the numbers of neutrino events (Rj? or R^f ) would contain 

1% or 0.2% error respectively, as -^p- or -^p- = ( 2 2 x x " i 6 ) — 1 • Experimentally 
obtained ratios of the numbers of the neutrino events above different energies from a 
large number of GRBs can be compared with their theoretically possible values Rj? or 
R^f . Detection of neutrinos from a large number of GRBs would reduce the statistical 
error in the measurement of the ratios of the numbers of neutrino events from them. 
If the systematic error of the detector in the measurement of number of neutrino 
events is 4% then the Dirac/ Majorana nature of neutrino would be distinguishable in 
that detector for the decay scenario with both i/ 2 and z/3 unstable using GRBs with 
Aa 2 ~ 0.001 as sources of ultrahigh neutrinos. 
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6 Conclusion 



We have calculated the ratios of the expected numbers of vertically upward going 
v T events from a single GRB in Multi-OWL experiment above energies 3 x 10 15 and 
2.2 x 10i m eV, including the oscillation and non-radiative two body decay of neutrinos. 
The ratios of the numbers of active neutrino events in the detector would be different 
depending on whether neutrino is a Dirac or a Majorana particle. These theoretically 
calculated ratios carry errors only from the error in the determination of CK2, which 
is obtained from the low energy photon data of the GRB. Measurement of the ratios 
of the numbers of neutrino events from a large number of GRBs would reduce the 
experimental statistical error. If the systematic error in detecting neutrino events is 
4%, and the error in determination of a<i is of the order of 0.001, then we expect the 
Dirac/Majorana nature of neutrinos to be revealed in Multi-OWL for the scenario of 
neutrino decay when both v<i and v$ are unstable. 
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APPENDIX A 

The parameters of a Gamma Ray Burst 



GRBs are short lived bursts of 7-rays of durations milliseconds to several minutes. 
After a burst emission of X-rays, optical and radio waves continue for several hours 
to months. This is called the afterglow phase of a GRB. The observed radiation is 
produced both during the burst phase and the afterglow phase. In the fireball model 
GRBs are produced by the dissipation of the kinetic energy of a relativistic expanding 
fireball. keV-MeV energy photons are expected to be produced by synchrotron shock 
accelerated electrons. A GRB occurs at a redshift "z", it has a Lorentz factor £ T" 
and it is expected to emit energy of the order of 10 51 erg in keV to MeV energy 7-ray 
emission. The observed photon spectrum on the Earth from a GRB can be expressed 
as a broken power law. The data obtained by BATSE (Burst and Transient Source 
Experiment) can be fitted using the parametrisation given by 



measured on Earth. The break energy is the characteristic observed energy of syn- 
chrotron photons produced by shock accelerated electrons. We denote the luminosity 
of a GRB in keV to MeV energy photon emission by "L 7 ". There is a detailed 
discussion on photon and neutrino production from GRBs in the fireball model in 
[26] . As the fireball expands relativistic shells come out from it with Lorentz factor 
T ~ 100 — 1000. The variability time of a GRB is the time interval between 
two successive emission of shells from it. When two shells collide with each other 
their kinetic energy is transformed into internal energy of the fireball. The fireball 
expands over a range of radii and reverse shocks are generated by the interaction of 
the expanding shells with the surrounding medium. Protons can be accelerated to 
very high energies by Fermi mechanism in internal and reverse shocks in the fireball. 
The low energy (keV-MeV) photon spectral break energy in the source rest frame is 
" e 7 ,Mev"> where e^ MeV = e^ obMeV (l + z), and it depends on the equipartition pa- 
rameters e e and e#, low energy photon luminosity of the GRB L 7 , Lorentz factor V 
and variability time t v . "e e " is the fraction of proton's shock thermal energy trans- 
formed to electrons and "e^" is the fraction of proton's shock thermal energy carried 
by the magnetic field in the fireball model. It is assumed that the wind luminosity 
carried by internal plasma energy L int , is related to the observed 7-ray luminosity 
through L int = £ 7 /e e . This assumption is justified because the electron synchrotron 
cooling time is short compared to the wind expansion time hence, electrons loose all 
their energy radiatively. The expression for photon spectral break energy can also be 



e 7j0 {, is the photon energy and e 7 
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expressed as 



rl/2 

1/2 3/2 -^7,52 



6 -y,MeV ~ 6 B e e -p 



2 



.t, 



2.5 t f,-2 

where, £ 7 ,52 = L 7 /10 52 , T 2 .5 = T/10 2 ' 5 and t„,-2 = t v /10~ 2 . Inside the fireball protons 
interact with photons to produce pions which subsequently decay to produce neutri- 
nos. There is a threshold energy for pion production in proton-photon interaction. 
For the energy of the proton-photon interaction to exceed the threshold energy for the 
A resonance, the proton energy in the source rest frame must satisfy the condition 

d, 

E b p > (1.4/< MeV )lt 5 x 10 7 GeV (9) 

The reference frame which has been defined as "observer's rest frame" in ^3], is 
referred as "source rest frame" in this paper. From the expression of proton spectrum 
break energy one can obtain the neutrino spectrum break energy. We assume that on 
the average 20% of the initial proton energy goes to the photo-pion. We also assume 
that the four final state leptons produced in the decay chain 7r + — > v^p^ — > v tl e + b> e v fl 
equally share the pion energy. Then the expression for neutrino spectrum break 
energy in the source rest frame is as follows 

r 2 

E b u = 7x 10 5 ^^GeV. (10) 

e -y,MeV 

The high energy pions produced due to proton-photon interactions inside the 
fireball are expected to radiate energy due to their synchrotron emission. One can 
compare the synchrotron energy loss time of a high energy pion with its lifetime for 
decay. Radiative energy loss of pions becomes important when their decay lifetime 
is larger than their synchrotron energy loss time. A break energy in the neutrino 
spectrum E s v would appear due to the synchrotron radiation emission of the pions 
before they decay to produce neutrinos. Assuming that a neutrino or an antineutrino 
takes only one fourth of the pion energy one can write down the expression for E s v 
from [TT] . 

El = 10 8 el/\ B 1/2 L-% 2 ri,t v> - 2 GeV (11) 

The total energy emitted by the GRB in electron, muon neutrinos and antineutrinos 
is E V ^ HE . 

rEvmax dN u 

E U) he = \ E u — r dE u (12) 

E l/m i n dUly 

where, u E um i n " and u E umax " are the minimum and maximum neutrino energy in the 
source rest frame. In the fireball model protons can be accelerated to a maximum 
energy of 10 21 eV. Since, a neutrino on the average receives 5% of the proton energy, 
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the maximum neutrino energy in the source rest frame u E umax " is approximately 
10 20 eV. 



13 



APPENDIX B 

Derivation of the factor W(E', E) 

We assume the decay of a heavy neutrino Vh of mass mh to a light neutrino v% of 
mass mi occurs through —> v\ + J, where J is a massless scalar (Majoron). The 
decay is given by the interaction Lagrangian density 

Li = g\V[v h J + g 2 Va$v h J + h.c. (13) 

The energy of the parent neutrino is E' and daughter neutrino is E. The total decay 
rate in the rest frame of the parent neutrino is 

2 2 

rLai = ~^fl\ Re ^\ 2 ( m i + m ^ + \^g 2 \ 2 (m h - mi f] (14) 

The above expression can be derived by standard methods [27J. In the labo- 
ratory frame the decay rate is T tota i(E') = ^T^ otal . For Majorana neutrinos, 
we have contribution from both conserved helicity and flipped helicity amplitudes 
Ftotai(E') = T f (E') + T nf (E'). The differential decay rate is dr totadE',E) _ The tQtal 
decay rate can also be written as 

W£ W E -^^V (15) 

2 

where, E m i n = —?E' and E max = E' . We work in the laboratory frame where 

h 

E,E' >> mh,mi. We obtain the energy distribution function of daughter neutrinos 
in the case of Majorana neutrinos as follows. 

rrr , rpt rp\ 1 dT tota l(E' , E) 

W M ^ m (E , E) = — W) . (16) 

The above expression can be expressed as a function of neutrino masses. 

2 2 

W Majorana (E>, E) = - mh with^E' <E<E>. (17) 
{mi - mf)E< mi 

For Dirac neutrinos only the daughter neutrinos, which conserve helicity, are con- 
tributing to the active species. So we define 

W Dirac (E', E) = 1 dT "f[ E '> E \ (is) 

and we obtain the energy distribution function as given below. 

W Dirac (E', E) = m h * J with—^-E < E < E . (19) 

\ m h ~ m i ) & m h 

Radiative decay going through transition magnetic moment will also lead to a 

similar result 28J. 
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APPENDIX C 

Derivation of the probability Prob{y T — > r) 

The energy of the u T is denoted by E Vt (0) before it enters the Earth. While v T is 
propagating through the Earth it is losing energy by neutral current interactions. The 
expression for energy loss of v T by neutral current interaction is 



(20) 



dx Jo ' ~'~ dl'. j (./•) 

= -p(x)N A a NC (E VT (x)){< E f (x) > -E Ut (x)) 

2X nc (E Ut (x)) 



We have used the fact that in each neutral current interaction the u T loses half of 
its energy [TE] in writing the final expression in eqn.(19). <tnc(E Vt (x)) is the cross 
section and \nc{E Vt (x)) is the mean free path of neutral current interaction. Ef(x) 
is the energy of the neutrino after it loses energy by neutral current interaction. The 
expression for the probability of neutral current interaction of v T is 

P^ UT (E Ur (0),E„ T (x),x)dE l/T (x) = 5(E Ur (x)-E Ur (0)e~ j ° ^cW)))dE UT (x) (21) 

The probability of a r production with energy between E T (x) and E T (x) + dE T (x) 
from a v T of energy E Vt (x) at a distance between x and x + dx is 



P^ T (E„M,E T (x),x)dE T (x)dx = p{x)N A dxdE T {x) d °cc{E» T {x) E T {x)) 

x P^ T (E I/T (0),E 1/T (x),x)dE UT (x). (22) 

where, <jcc{E Vt {x)) is the cross-section of charged current interaction for v T with 
energy E Ut (x), p(x) is the density of the medium at x and Na is Avogadro's number. 
We assume that the r lepton produced by charged current interaction carries 75% of 
the v T energy ^7]. The r leptons are losing energy during their propagation. There 
is discussion on energy loss of r leptons in j^I]. The energy loss by r leptons due to 
ionization is negligible compared to the other processes e.g. bremsstrahlung, e + e~ 
pair production and photonuclear processes for energies above 10 5 GeV. We have used 
the energy loss co-efficient of r lepton 

(3{E T {x)) = (1.6 + 6(E T {x)/10 9 GeV) - 2 ) x 10~ 7 g^cm 2 (23) 

The r lepton may decay or undergo charged current interaction to produce a v T . The 
v T produced from r decay carries 40% of the r lepton energy. The expression for the 
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probability of neutral current interaction including r decay and r charged current 
interaction is 



P„ r ^ T (E Vr (0),E UT {x),x)dE UT (x) = P^ 1/T (E UT (0),E UT (x),x)dE UT (x) 

+ dE VT {x) £ P^iE^xJ, E UT (x),x - xi) 

,dT T (E UT (x 1 ),E T (x 1 )) 



x 



dE VT { x \) 



drt ?^" T x (E u (x 1 ),EJx 1 ) , , Ar . (rn , ., 

+ ) % v V V 1JJ P (x 1 )N A )dE UT (x 1 )dx 1 

dE VT {xi) 



r [ " TX2 P Vr ^ T (E VT (0),E T (x 2 ),x 2 )dE T (x 2 )dx 2 
Jo Jo 

(24) 



The energy of the r lepton at x% and x 2 are related by the energy loss equation 
of r leptons. 

= -/3(£ T (x))£ T (a;) (25) 

ax 

The expression for the probability of charged current interaction is 

P Ut ^ t (E Ut (0), E t (x), x)dE T (x)dx = p(x)NAdE T (x)dx 

x / P Vt ^ t {E Ut (0),E Ut (x),x) 



dE T (x) 

The r leptons produced by charged current interaction lose energy before emerg- 
ing from the surface of the Earth. The range of r is approximately R T = 
1/ j3(E T (x)) \n(E T (x) / E t h) from eqn.(24). E t h is the threshold energy of the detec- 
tor. The range of tau leptons can be written as a function of E Ut (x) by replacing 
E T (x) with 0.75£ , j /t (x). The final expression for the probability of emission of r lep- 
tons from the surface of the Earth having energies above the threshold energy of the 
detector is obtained after integrating eqn.(25) over E T (x) and x. 

Prob(u T r) = P Vt ^ t {E Ut {0), 0.75E VT (x),x)dx (27) 
Jo 

We have solved eqn.(23) and (25) iteratively to find the value of the probability 
on the left hand side of eqn.(25). We use the charged current and neutral current 
interaction cross-sections from [2H] and the Earth's density profile from [3U]. Our 
result converges after two iterations. Finally we have use the probability calculated 
in eqn.(26) Prob[y T — ► r) for our calculation of number of r events in the detector. 
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